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Abstract 

Background and Purpose: Little is known about the interactions between the default mode network (DMN) subregions in 
relapsing-remitting multiple sclerosis (RRMS). This study used diffusion tensor imaging (DTI) and resting-state functional MRI 
(rs-fMRI) to examine alterations of long white matter tracts in paired DMN subregions and their functional connectivity in 
RRMS patients. 

Methods: Twenty-four RRMS patients and 24 healthy subjects participated in this study. The fiber connections derived from 
DTI tractography and the temporal correlation coefficient derived from rs-fMRI were combined to examine the inter- 
subregion structural-functional connectivity (SC-FC) within the DMN and its correlations with clinical markers. 

Results: Compared with healthy subjects, the RRMS patients showed the following: 1) significantly decreased SC and 
increased FC in the pair-wise subregions; 2) two significant correlations in SC-FC coupling patterns, including the positive 
correlation between slightly increased FC value and long white matter tract damage in the PCC/PCUN-MPFC connection, 
and the negative correlations between significantly increased FC values and long white matter tract damage in the PCC/ 
PCUN-bilateral mTL connections; 3) SC alterations [log(N track) of the PCC/PCUN-left IPL, RD value of the MPFC-left IPL, FA 
value of the PCC/PCUN-left mTL connections] correlated with EDSS, increases in the RD value of MPFC-left IPL connection 
was positively correlated to the MFIS; and decreases in the FA value of PCC/PCUN-right IPL connection was negatively 
correlated with the PASAT; 4) decreased SC (FA value of the MPFC-left IPL, track volume of the PCC/PCUN-MPFC, and log(N 
track) of PCC/PCUN-left mTL connections) was positively correlated with brain atrophy. 

Conclusions: In the connections of paired DMN subregions, we observed decreased SC and increased FC in RRMS patients. 
The relationship between MS-related structural abnormalities and clinical markers suggests that the disruption of this long- 
distance "inter-subregion" connectivity (white matter) may significantly impact the integrity of the network's function. 
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Introduction 

Pathological [1-4] and neuroimaging [5—7] studies have 
demonstrated the widespread involvement of the gray matter 
and white matter of the brain in multiple sclerosis (MS). In MS, 
minimal initial structural damage can be sufficient to trigger a 
strong response in function and hyperactivity [8,9]. Structural 
damage limits functional reorganization [10] until burnout, 
resulting in hypoactivity and dysfunction [11]. Although the 
structures has been considered as the underlying foundation of 
functional change. But in fact, it is controversial in the relationship 
between the cerebral structure and function or clinical function 
[12—14]. Understanding the pathophysiological mechanisms 
relating structural damage and functional alterations in MS and 



its correlations with clinical indices remains a challenge in clinical 
research [15-17]. 

Recent studies of MS have reported significant synchronization 
or connectivity alterations in diverse disease stages [8,9,18,19], 
including changes in the default-mode network (DMN) [11,18— 
20], sensorimotor network (SMN) [21,22] and frontoparietal 
network (FPN) [9,23] (reviewed in [24]). The DMN, a recent focus 
of functional MRI studies of MS, may be particularly important 
for understanding the impact of the disease on the brain [25,26]. 
Especially in relapsing-remitting MS (RRMS), different alterations 
within the DMN have been described, including decreased [27], 
increased [19] and mixed [18] changes as well as no significant 
differences [8]. Widespread functional abnormalities between the 
DMN and other large-scale networks [28] have also been 
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reported, but the correlation between significandy altered brain 
activity and clinical characteristics/neuropsychological tests scores 
is inconsistent [24]. 

The DMN comprises multiple dissociated subregions that are 
engaged to facilitate functional construction related to cognition 
through the long-distance connectivity "backbone". In contrast to 
local interactions, this long-distance cortico-cortical (inter-subre- 
gion) interaction is involved in the global communication of the 
whole network [29]. The anatomical and functional connectivity 
(FC) linking to paired DMN subregions has been accessed using 
DTI tractography and resting-state functional connectivity 
[13,14,30-33]. 

However, little is known about the connectivity "backbone" 
linking the diverse DMN subregions in MS. Studying the long- 
distance ("backbone") inter-subregion connections may provide 
more information that is critical to fully understanding MS-related 
functional alterations. In this study, we hypothesized that both 
structural connectivity (SC) and FC in paired DMN subregions are 
disrupted, as are its correlations with clinical indices in RRMS. To 
test this hypothesis, we combined probabilistic tractography and 
pair-wise FC to characterize the connectivity of long white matter 
tracts and the synchrony of intrinsic neuronal activity in paired 
DMN subregions and to examine the mechanisms that might 
underlie the observed alterations. Studying the relationship 
between the SC and FC of anatomically separate DMN 
subregions may enrich our understanding of the neural under- 
pinnings of RRMS. 

Materials and Methods 

Participants 

We recruited 24 patients with clinically definite MS at the First 
Affiliated Hospital of Nanchang University from May 2010 to 
December 2013, according to McDonald's criteria [34]. The 
inclusion criteria for the patients were as follows: RRMS course 
[35]; Expanded Disability Status Scale (EDSS) score <2.5 
(corresponding to minimally disabled [36]); and treatment with 
immunomodulatory medication (20 with fj-interferons, 4 with 
Glatiramer acetate). None of the recruited patients had any 
relapses or cortico-steroid treatment during the month preceding 
the MR acquisition. Twenty-four healthy control participants 
from the local community were matched to the patients 
individually for sex, age, and education level. Their clinical and 
demographic data are shown in Table 1. The present study was 
approved by the Medical Research Ethics Committee and the 
Institutional Review Board of the First Affiliated Hospital of 
Nanchang University. All subjects signed written consent forms for 
participation in the study. 

Image acquisition 

All participants were scanned using a 3.0 T MRI scanner (Trio 
Tim, Siemens Medical Systems, Erlangen, Germany). T 2 W, T\W, 
DTI and resting-state fMRI (rs-fMRI) images were acquired using 
the following sequences: 1) T 2 -weighted turbo spin echo imaging 
(repetition time [TR] /echo time [TE] = 5 1 00/ 1 1 7 ms; number of 
excitations [NEX] = 3; echo train length =11; matrix = 416x416; 
field of view [FOV] = 240 x240 mm; slice number = 22; slice 
thickness = 6.5 mm; and orientation = axial); 2) T\ -weighted 3-D 
imaging (TR/TE = 1900/2.26 ms; NEX= 1; matrix = 240 x 256; 
FOV = 215x230 mm; slices number = 176; slice thick- 
ness = 1 .0 mm; and orientation = sagittal); 3) spin echo single-shot 
echo planar imaging (TR/TE = 7200/104 ms; NEX = 2; ma- 
trix = 128 x 128; FOV = 230 x230 mm; slices number = 49; slice 
thickness = 2.5 mm; orientation = axial; 64 nonlinear diffusion 



weighting gradient directions with b = 1000 s/mm 2 and 1 addi- 
tional image without diffusion weighting [i.e., b = 0 s/mm 2 ]); and 
4) rs-fMRI scan using an echo planar imaging (EPI) sequence with 
the foUowing parameters: TR/TE = 2000/30 ms; flip angle = 90°; 
FOV = 200x200 mm; matrix = 64x64; 30 interleaved axial slices 
with 4-mm thickness with an interslice gap of 1.2 mm; and 
number of time points = 240. During rs-fMRI scanning, subjects 
were instructed to keep their eyes closed, not to think about 
anything in particular, and not to fall asleep. A foam pad was used 
to minimize the head motion of all subjects. 

Functional MRI data preprocessing 

The rs-fMRI data were preprocessed using the SPM8 software 
package (http://www.fil.ion.ucl.ac.uk) and the Data Processing 
Assistant for Resting-State fMRI Advanced Edition (DPARSFA) 
V2.2 (http://www.restfmri.net) running in Matlab 2012a (Math- 
works, Natick, MA, USA). This process involved removing the first 
1 0 time points, slice timing, voxel-specific head motion calculation 
and correction to adjust the time series of the images, using unified 
segmentation of the Ty-weighted images for non-linear registration 
and normalization of the functional images to the Montreal 
Neurological Institute (MNI) space with 3x3x3 mm 3 re-sampling, 
then spatial smoothing was performed using a 6-mm full-width- 
half-maximum Gaussian kernel. We used temporal band-pass 
filtering (0.01 Hz<f<0.08 Hz) and multiple regression analysis of 
nuisance variables from the BOLD data, which included a 
ventricular signal averaged from ventricular regions of interest 
(ROIs), a white matter signal averaged from white matter ROIs, a 
whole brain signal averaged across the whole brain, six head 
realignment parameters obtained by rigid body head motion 
correction, the derivatives of each of these signals and a voxel- 
specific head motion parameter [39]. 

Default mode network subregion extraction 

The default mode network subregions were extracted using a 
similar method to that described in a previous study [14], which is 
briefly described as follows: 

1. Group independent component analysis (ICA) was first 
performed using the Group ICA of fMRI Toolbox (GIFT) 
(http://icatb.sourceforge.net/, vision 2.0e) with the Infomax 
algorithm, which decomposed the smoothed data of each 
individual in the RRMS and healthy control groups into 63 
and 69 independent components (estimated by the minimum 
description length (MDL) criterion [40]), respectively. The 
RRMS and healthy control groups in this study were estimated 
and analyzed separately to prevent the mixing of the specific 
resting-state network pattern from each group [41]. 

2. The largest spatial correlation with DMN templates was used 
to select the components to be retained for further analysis. 
The DMN templates form the Medical Image Analysis (MIA) 
laboratory (http://mialab.mrn.org/index.html), which was 
divided into six subregions: bilateral inferior parietal lobules 
(IPLs), bilateral medial temporal lobes (mTLs), posterior 
cingulate cortex/precuneus (PCC/PCUN) and medial pre- 
frontal cortex (MPFC). Each of the subregions was manually 
divided into ROIs based on the union of the DMN maps from 
the RRMS and healthy control groups (REST VI. 8, http:// 
www.restfmri.net). Detailed information on ROIs is shown in 
Figure 1 and Table S 1 . 

3. These six subregions were used for subsequent SC and FC 
analyses in the patient and control groups, respectively. 
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Table 1. Demographics and clinical characteristics of the control subjects and RRMS patients. 







RRMS patients (n = 24) 


Control subject (n = 24) 


P values 


Gender (M/F) 


8/16 


8/16 


>0.99 


Mean age (range) (years) 


39.5 (20-56) 


39.6 (21-56) 


0.99 


Mean disease duration (range) (months) 


34.1 (3-187) 






Mean TWMLL* (range) (ml) (normalized) 


18.17 (0.43-79.41) 






Mean BPF (range) 


0.828 (0.78-0.86) 


0.853(0.82-0.88) 


0.000 


Median EDSS (range) 


1.604 (1.0-2.5) 


0 




Mean PASAT (range) 


84.0 (61-103) 


97.87 (82-118) 


0.002 


Mean MFIS-5 (range) 


11.29 (6-17) 


0.29 (0-1) 


0.000 



Note: M = male; F = female; TWMLL = total white matter lesion loads; BPF = brain parenchymal fraction; EDSS = expanded disability status scale; PASAT = paced auditory 
serial addition test; MFIS = modified fatigue impact scale; 

*The measurement procedures for TWMLL and BPF in the RRMS patients have been previously described by [37] and [38], respectively (see Figure SI and Supporting 
Information SI). 

doi:1 0.1 371 /journal.pone.01 01 1 98.t001 



Calculation of synchrony in the intrinsic neuronal activity 
of pair-wise subregions 

To extract and compute the time series of each pair of 
subregions from two groups, we used a similar method to that 
described in previous studies [32,41], including multiple regression 
of nuisance variables and band-pass filtering (0.01-0.08 Hz). 
Temporal correlation coefficients were calculated only between 
each pair of regions using the Pearson correlation to quantify the 
FC values (Figure 2). 

Characterizing the connectivity of long white matter 
tract from DTI data 

To investigate the SC of pair-wise subregions, the group-space 
ROI was coregistered to individually (native) DTI space using 
FLIRT methods [42] . They are briefly described as follows: 

1. Using FMRIB's diffusion toolbox (FDT v2.0, http://fsl.fmrib. 
ox.ac.uk/fsl/fslwiki/FDT), DTIFIT was used to fit a single 
tensor model at each voxel of the preprocessed eddy current 
corrected diffusion weighted data. 

2. Probabilistic tractography was performed using FMRIB's 
diffusion toolbox (FDT v2.0, http://fsl.fmrib.ox.ac.uk/fsl/ 
fslwiH/FDT). BEDPOSTX was used to model 5000 iterations 
within each voxel with a curvature threshold of 0.2, a step 
length of 0.5 and a maximum number of 2000 steps [43]. Pair- 




Figure 1. The union of the DMN mask identified and extracted 
with independent component analysis from RRMS and healthy 
control groups (LH: left hemisphere; RH: right hemisphere). 

doi:1 0.1 371 /journal.pone.01 01 1 98.g001 




FC analysis SC analysis 



Figure 2. A flowchart of calculating the correlation between 
structural connectivity and functional connectivity within the 
default mode network (DMN) subregions. (1) DMN extraction was 
performed by ICA. (2) The subregions were extracted as regions of 
interest (ROIs) within the DMN. (3) The time series in each ROI was 
extracted with the DMN. (4) The temporal correlation coefficients 
between each ROI pair were quantified within the DMN. (5) 
Preprocessed motion and eddy current distortion correction of the 
model distributions of the relevant parameter (Monte Carlo sampling) 
was performed using FMRIB's diffusion toolbox (FDT v2.0). (6) Pair-wise 
ROI (AND logic) probabilistic tractography was used to calculated the 
distribution of fiber orientations. (7) Features of the long white matter 
tracts (above threshold 0.2) remaining fiber bundles connecting each 
pair of ROIs were compared between the two groups, including the 
strength of a pathway and the above-threshold standard DTI 
parameters. 

doi:1 0.1 371/journal.pone.01 01 1 98.g002 
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wise ROIs ("AND" logic) were used to calculate a distribution 
of fiber orientations. 

3. The connectivity of long white matter tracts was at a 
normalized probability value of 0.20 [31,44] and visually 
inspected to confirm successful tracing in each individual by 
multiple experienced technicians. 

4. Features of the long white matter tracts (above threshold) 
connecting each pair of ROIs were compared between the two 
groups. The features included the strength of a pathway, 
indicated by the volumes and mean track count per WM 
region of the paired ROIs [43], and the above threshold 
standard DTI parameters, including fractional anisotropy (FA) 
and three diffusivity measurements [average mean diffusivity 
(MD), axial diffusivity (AD) and radial diffusivity (RD)] . 

Comparison of structure and function 

We assessed the character of the SC and FC measures in three 
steps: 1) structural and functional connections were identified of 
pair-wise DMN subregions, specifically, the fiber connections as 
being either present or absent, identified by visual inspection; 2) 
the connectivity strength of pair-wise subregions was determined 
and inter-groups compared separately for FC and SC; and 3) the 
relationship of SC-FC coupling was determined by linear 
correlation analysis for each pair-wise connection. Pair-wise 
probabilistic tractography tracks in which a significant difference 
was detected between groups were excluded from this analysis. 

Statistical analysis 

For the FC and SC separately, we compared the 15 pairs FC 
between the ROIs (i.e., PCC/PCUN-MPFC, MPFC-right IPL) 
and eight pairs SC using repeated one-way repeated-measures 
analysis of variance (ANOVA) (SPSS Inc., Chicago, IL, USA). 
Post-hoc T-tests were used to identify differences in SC and FC 
between two groups (P<0.05, Bonferroni corrected). Linear 
regression analysis was performed using SPSS 13.0 (SPSS Inc., 
Chicago, IL, USA) to investigate the relationship between 
neuropsychological tests scores and damage to SC/FC and the 
effect on clinical parameters and indices within RRMS patients. 
The threshold was set at a significance level of P<0.05 and was 
corrected for multiple comparisons using the Bonferroni correc- 
tion. 

Results 

Default mode network subregion extraction by ICA 
analysis 

In the present study, DMN maps revealed a typical spatial 
pattern in both the patient and control groups (P<0.05, FDR 
corrected), and six subregions were extracted from the union of the 
DMN maps of two groups (Figure 1). The details of the brain 
regions in the DMN subregions are reported in the supporting 
materials (Table SI). 

Functional connectivity measures of pair-wise subregion 
coherence 

Seed-based FC was then used to measure the FC strength for 
each pair-wise connection between the DMN subregions. Figure 3a 
presents the mean FC coefficients of the pair-wise subregions. An 
analysis of the FC coefficients of seven paired subregions (MPFC- 
left IPL, PCC/PCUN-left mTL, PCC/PCUN-right mTL, left 
IPL-left mTL, left IPL-right mTL, right IPL-right mTL, and left 




Figure 3. A set of matrices of functional and structural analyses 
of pair-wise default mode network subregions in the RRMS 
patient and healthy control (HC) groups. The figure (a) is a 
correlation matrix of average time series of subregions (ROIs) in the 
network (subregion labels are on the left and bottom). The figures (b-s) 
show the strength of the structural connectivity matrix [of average log 
(N track) and volumes] of each pair-wise connection. The figures (d-g) 
are above-threshold (0.2) standard DTI parameters: including the 
fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD) 
and radial diffusivity (RD) of each pair-wise connection. The black arrow 
indicates increased and the white arrow decreased alteration in the 
patients with RRMS. Notably, structural connectivity was only compared 
without the dashed boxes. 
doi:1 0.1 371/journal.pone.01 01 1 98.g003 

mTL-right mTL) showed significant differences between the two 
groups at P<0.05 (ANOVA; Tables S2 and S3). 

Structural connectivity measures of long white matter 
tracts between pair-wise subregions 

Figure 4 showes the long white matter tracks identified by 
probabilistic tracking between pair-wise subregions. Compared 
with the healthy control group, a lower probabilistic track number 
was detected in the RRMS group (Table S3 and examples in 
Figure S2). 

In this study, we only compared the structural connectivity 
measures of non-significant differences in detectable tractography 
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Right Pos. 

HC 



Left Ant. 
Right Pos. 



between the PCC/PCUN and the MPFC 



RRMS 



Figure 4. An example of the long white matter fibers of pair- 
wise subregions within the default mode network detected by 
probabilistic tractography in one control subject (a) and one 
RRMS patient (c) in the native space. Figures (b) and (d) 
correspond to 3-dimensional probabilistic tractography. (red: fiber 
tract, blue: mask of subregions, Sup = Superior, Inf = Inferior, Ant = An- 
terior, Pos = Posterior). 
doi:1 0.1 371 /journal.pone.01 01 1 98.g004 



(eight connections) in the pair-wise subregions between two 
groups. The SC strength was assessed by measuring the volumes 
and mean track count along the above-threshold reconstructed 
tracks (Figures 3b, 3c and Tables S2, S3). In the patients with 
RRMS, smaller volumes were detected in the track linking the 
PCC/PCUN and the MPFC (4.55 xlO 3 mm 3 vs. 2.26 xlO 3 mm 3 ) 
and the track linking the MPFC and the right IPL (3.88 x 1 0 3 mm 3 
vs. 2.17 xlO 3 mm 3 ); lower counts [Log(N track)] were detected in 
the track linking the PCC/PCUN and the MPFC (5.69 vs. 5.21), 
the track linking the PCC/PCUN and the left mTL (4.77 vs. 2.83), 
and the track linking the PCC/PCUN and the right mTL (4.90 vs. 
2.84). 

The integrity of the above-threshold reconstructed tracks was 
assessed by measuring the standard DTI parameters (Figures 3d- 
3g and Tables S2, S3). In the patients with RRMS, decreased FA 
and increased diffusivity measurement (MD, AD and RD) values 
were detected in all eight tracts (P<0.05, Table S3), with the 
exception of the AD values of the track linking the PCC/PCUN 
and the bilateral IPL, and the track linking between the bilateral 
IPLs. 

Relationship between structure and function measures 

Linear regression analyses (Figure 5) revealed that in RRMS 
patients, increases in the MD and AD values of the track linking 
the PCC/PCUN and the MPFC were positively correlated with 
slight increases in FC values (Figure 5a, r= 0.545, P— 0.006; 
Figure 5b, r— 0.446, P— 0.029, respectively). The smaller decrease 
in the volumes of the track linking the PCC/PCUN and the left 
mTL was positively correlated with increased FC values (Figure 5c, 
r— 0.517, P— 0.014). An increase in the MD values of the track 
linking the PCC/PCUN and the right mTL was negatively 
correlated with an increase in the FC value (Figure 5d, r= —0.453, 
P— 0.034). Figure S3 shows an example of the relationship 
between the FC and SC measures in others threshold (0.05, 0.10, 
0.15, 0.20, 0.25, and 0.30) of probabilistic tractography. 
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Figure 5. The Pearson correlation between structural connec- 
tivity (SC) indices (y axis) and functional connectivity (FC) 
indices (x axis) of pair-wise subregions within the DMN of 
RRMS patients. The FC value positively correlated with the MD value 
(a) and FC value (b) of the tract between the PCC/PCUN and the MPFC 
in RRMS patients (P<0.05, Bonferroni corrected). The FC values was 
positively correlated with the volume (c) of the tract between the PCC/ 
PCUN and the left mTL in RRMS patients (P<0.05, Bonferroni corrected). 
The FC value positively correlated with the MD value (d) of the tract 
between the PCC/PCUN and right mTL (P<0.05, Bonferroni corrected). 
AD = axial diffusivity, IPL = inferior parietal lobule, MD = mean diffusivity, 
mTL = medial temporal lobes, MPFC = medial prefrontal cortex, PA- 
SAT = paced auditory serial addition test, PCC/PCUN = posterior cingu- 
late cortex/precuneus, RRMS = relapsing-remitting multiple sclerosis. 
doi:1 0.1 371/journal.pone.01 01 1 98.g005 

Linear regression results between abnormal connectivity 
indices and clinical markers of RRMS 

Linear regression analyses revealed that in RRMS patients, 
increases in FC values between the PCC/PCUN and left mTL 
were negatively correlated with the PASAT score (Figure 6a, 
r= -0.435, P= 0.033). 

The mean track count of the track linking the PCC/PCUN and 
the left IPL was negatively correlated with the EDSS value 
(Figure 6b, r= -0.420 P= 0.041). A decrease in the FA value of 
the track linking the PCC/PCUN and the left mTL was negatively 
correlated with the EDSS value (Figure 6c, r- -0.434, P= 0.044). 
An increase in the RD value of the track linking the MPFC and 
the left IPL was positively correlated with the EDSS and MFIS 
values (Figure 6d, r= 0.539 P= 0.007; Figure 6e, r= 0.410, 
P— 0.046, respectively). A decrease in the FA value of the track 
linking the PCC/PCUN and the right IPL was positively 
correlated with the PASAT value (Figure 6f, r= 0.447, P= 0.028). 

A decrease in the FA value of the track linking the MPFC and 
the left IPL was positively correlated with the BPF value (Figure 6g, 
r— 0.430 P= 0.036). A decreased volume of the track linking the 
PCC/PCUN and the MPFC was positively correlated with the 
BPF value (Figure 6h, r- 0.473 P= 0.020). The mean track count 
of the track linking the PCC/PCUN and the left mTL was 
positively correlated with the BPF value (Figure 6i, r— 0.478 
P= 0.024). 

Discussion 

In the present study, we showed a significandy decreased inter- 
subregion SC in the DMN and increased FC of pair-wise 
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Figure 6. The Pearson correlation between structural/functional connectivity indices (x axis) and clinical markers (y axis) of pair- 
wise subregions within the DMN of RRMS patients. The functional connectivity (FC) value (a) between the PCC/PCUN and the right mTL 
negatively correlated with the PASAT in RRMS patients (P<0.05, Bonferroni corrected). The EDSS score negatively correlated with the log (N track) of 
the tract between the PCC/PCUN and the left IPL (b) and negatively correlated with the FA value of the tract between the PCC/PCUN and the left mTL 
(c) in RRMS patients (P<0.05, Bonferroni corrected). The EDSS score positively correlated with the MD value of the tract between the MPFC and the 
left IPL (d) in RRMS patients (P<0.05, Bonferroni corrected). The PASAT score positively correlated with the FA value of the tract between the PCC/ 
PCUN and the right IPL (e) in RRMS patients (P<0.05, Bonferroni corrected). The MFIS score negatively correlated with the MD values of the tract 
between the MPFC and the left IPL (f) in RRMS patients (P<0.05, Bonferroni corrected). The BPF score positively correlated with the FA value of the 
tract between the MPFC and the left IPL (g) in RRMS patients (P<0.05, Bonferroni corrected). The BPF score positively correlated with the volume of 
the tract between the PCC/PCUN and the MPFC (h) in RRMS patients (P<0.05 : Bonferroni corrected). The BPF score positively correlated with the log 
(N track) value between the PCC/PCUN and the left mTL (i) in RRMS patients (P<0.05, Bonferroni corrected). AD = axial diffusivity, EDSS = expanded 
disability status scale, FA = fractional anisotropy, IPL = inferior parietal lobule, MD = mean diffusivity, mTL = medial temporal lobes, MPFC = medial 
prefrontal cortex, PASAT = paced auditory serial addition test, PCC/PCUN = posterior cingulate cortex/precuneus, RRMS = relapsing-remitting multiple 
sclerosis, TWMLL = total white matter lesion loads. 
doi:1 0.1 371 /journal.pone.01 01 1 98.g006 



subregions in minimally disabled RRMS patients. We also found 
two significant correlations between the FC and SC measures: 
slightly increased FC value was positively correlated with long 
white matter tract damage in the PCC/PCUN-MPFC connection, 
and significantly increased FC values were negatively correlated 
with long white matter tract damage in the PCC/PCUN-bilateral 
mTL connections. Relationships between MS-related structural 
abnormalities and clinical markers were detected more often than 
relationships between functional abnormalities and clinical mark- 
ers. This study, linking the increased inter-subregion functional 
coherence and decreased long white matter connections of DMN 
subregions, provides new insights into the functional reorganiza- 
tion in RRMS. 



Decreased connectivity of long white matter tracts 
between the default mode network subregions 

DTI measures (region of interest analysis [6,45] and tract-based 
spatial statistics [46], for example) are sensitive to white matter 
damage in MS, including interior focal lesions and normal 
appearing white matter, and these regional DTI alterations are 
related to disability and cognitive deficits [47]. Although damage 
to the integrity of fiber tracts has been investigated in MS, very 
littie is known about alterations in the fiber connections between 
the distant subregions within the DMN. The long white matter 
tracts between subregions are the structural highways of the DMN, 
and they enable information to travel quickly from one subregion 
to another subregion to facilitate the ongoing interregional 
neuronal communication [14]. In contrast to deterministic 
tractography (just a single trajectory reconstruction), probabilistic 
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tractography was performed in our study to investigate the 
strength of SC and the above-threshold standard DTI parameters 
in the long white matter tracts. Probabilistic tractography methods 
are used to propagate a large number of pathways passing through 
the seed point, and the pathway orientations are drawn from a 
distribution of possible orientations [43,44]. The probabilistic 
tractography methods can be used to (facilitate) easily reconstruct 
the existence of direct neuroanatomies] connections between 
DMN subregions [14,30]. 

The potential of probabilistic tractography to detect abnormal- 
ities in the long white matter tracts of patients with MS has been 
suggested in recent studies [48-51]. These studies have reported 
the sensitive detection of tissue damage associated with cervical 
cord relapse [48] and damage to the corticospinal tract [49-51] 
and the transcallosal hand motor fibers [49], suggesting axonal 
degeneration and myelin breakdown in probabilistic white matter 
tracts, these measures correlate with disability [50,51]. In our 
study, probabilistic tractography analysis showed decreases in SC 
strength and structural damage derived from standard DTI 
parameters (decreased FA and increases in three diffusivity 
measurements) in the tracts connecting the DMN subregions. 
These findings further support the existence of inter-subregion 
anatomical disconnection and microstructural damage (axonal 
and myelin loss) of entire long white matter tracts in the DMN 
[19]. 

Notably, with significant anatomical disconnection (without 
white matter connection) in the DMN of patients with RRMS, the 
SC strength and the derived standard DTI parameters were 
statistically compared only for connections within eight paired 
subregions showing no significant indifferences. 

Increased functional connectivity in the pair-wise default 
mode network subregions 

In the present study, increased functional connectivity was 
observed in the patients with RRMS using pair-wise subregion 
coherence analysis, although not for all paired connections. 
Increased functional connectivity was detected in seven paired 
subregions (MPFC-left IPL, PCC/PCUN-left mTL, PCC/PCUN- 
right mTL, left IPL-left mTL, left IPL-right mTL, right IPL-right 
mTL, and left mTL-right mTL) within the DMN. Pair-wise base 
coherence analysis provides a simply and/ or straightforward FC 
metric [52]. Recently, various MS-related FC alteration patterns 
have been reported using rs-fMRI, including decreased [27], 
increased [19] and mixed [18] changes as well as no significant 
differences [8] . Our results may suggest an increased neural effort 
in cortical activity due to the presence of structural damage the 
connecting pathways in minimally disabled RRMS patients [19]. 

Relationship of SC-FC coupling measures 

Functional reorganization might be an important pathophysi- 
ological factor in controlling information processing deficiencies in 
relatively early MS. Neuronal fiber tracts are bundles that have an 
extremely large number of axons that act as a bridge connecting 
two larger subregions over long distances [53], known as the 
connectivity "backbone". In general, inter-subregions that are 
clearly and directly linked by anatomical long white matter tracts 
should exhibit a high level of functional communication [13]. The 
relationship pattern has previously investigated the use of both SC 
and FC analyses, which have suggested that the functionally linked 
DMN subregions reflect the underlying structural connectivity 
architecture [13,14]. These findings have enriched our under- 
standing of the neuro-pathophysiological mechanisms of disease 
[32,41]. 



In this study, significant relationships between SC-FC coupling 
measures were observed in the connection among the PCC/ 
PCUN and the MPFC (MD values vs. FC coefficients, r= 0.545, 
P= 0.006; AD values vs. FC coefficients, r= 0.446, P= 0.029, 
respectively), the connection between the PCC /PCUN and left 
mTL (volumes of the fiber tract vs. FC coefficients, r= 0.571, 
P= 0.014) and the connection between the PCC/PCUN and right 
mTL (MD values vs. FC coefficients, r= -0.453, P= 0.034). We 
identified three different MS-related relationships in SC-FC 
coupling measures: a) a slightly increased FC was positively 
correlated with damage to SC; b) significantly increased FC was 
negatively correlated with damage to SC; and c) dissociation of 
SC-FC coupling measures. The probabilistic tractography and 
resting-state functional connectivity approach is advantageous for 
characterizing the relationship of SC-FC coupling [31,54]. In 
minimally disabled MS patients, increased synchronization in the 
DMN have been observed [8], revealing a compensatory 
mechanism in response to structural damage (relationship pattern 
a), which seems to be finite. At some point, structural damage will 
limit the functional compensation [19] (relationship pattern b), but 
the increase in MS-related SC-FC coupling falls somewhere 
between these two patterns (relationship pattern c). In the 
dissociation of SC-FC coupling pattern, a normal positive 
correlation between the SC and FC measures does not exist 
[3 1 ,54] , even a tendency toward a negative correlation (Figure S3). 
However, to our knowledge, no studies have attempted to link the 
SC and FC of inter-subregions within the DMN in RRMS patients 
before this study. 

Relationship between abnormal connectivity indices of 
inter-subregions and clinical markers of RRMS 

This study provided strong evidence in support of a structural 
disconnection in patients with RRMS. Our study highlighted that 
altered SC in the long white matter tract [e.g., log (N track) of the 
PCC/PCUN-left IPL connection, RD values of the MPFC-left 
IPL connection, FA values of the PCC/PCUN-left mTL 
connection] correlated with the EDSS, that increases in the RD 
value of the MPFC-left IPL connection positively correlated with 
the MFIS and that decreases in the FA value of the PCC/PCUN- 
right IPL connection negatively correlated with the PASAT. The 
EDSS is used to quantify and monitor disability changes in MS. A 
decreased PASAT indicates slowed cognitive processing. The 
fatigue-related MFIS-5 scale is an abbreviated version used to 
assess the impact of fatigue on physical, cognitive, and psychoso- 
cial functioning. Previous studies have demonstrated that fatigue, 
EDSS severity and lower PASAT scores in RRMS patient are 
correlated with the white matter integrity [55,56] . In this study, we 
determined that structural damage to the inter-subregions within 
DMN was associated with clinical markers. 

Another highlight of our study was that SC abnormalities of the 
inter-subregions within the DMN in RRMS patients were 
significantly correlated with brain atrophy, rather than the total 
white matter lesion load. Lesion volume is considered to have little 
correlation with clinical scales when compared with atrophy [38], 
and our result conformed this opinion. 

Furthermore, the FC value of the PCC/PCUN-right mTL was 
negatively correlated with the PASAT score. Previous studies have 
demonstrated that the clinical severity of RRMS is inversely 
correlated with functional connectivity within the DMN [9,23]. In 
this study, increases in the FC value of the PCC/PCUN-right 
mTL were also negatively associated with myelin loss (higher MD 
values). This association likely occurred because structural damage 
limits the functional compensation, although functional reorgani- 
zation was ineffective in the PCC/PCUN-right mTL connection. 
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Study Limitations 

Importantly, large ROIs could blur the results of functional 
connectivity analyses. In this study, we denned ROIs that comprised 
more than one brain area based only on the union of the DMN maps 
from the RRMS and healthy control groups and without the use of 
anatomical information or prior seed coordinates from the literature. 
The selection of the DMN subregions remains controversial in SC-FC 
coupling studies. Another limitation to this study could be that the DTI 
probabilistic tractography analyses were restricted to the pair-wise 
connections. Additionally, there are also a few limitations to probabilistic 
tractography including the possibly highly reproducible topology of the 
reconstructed pathways [31,43] and lack of a statistical consensus on the 
probabilistic tractography threshold [44]. Moreover, it was demon- 
strated that robust FC was present between regions not linked by 
cortico-cortical fiber projections. Limited assessment of cognitive 
function was implemented in this study. Finally, an explorative study 
of the altered SC and FC of inter-subregions is warranted in the future 
to confirm or supplement the findings of this study. 

Conclusions 

In the present study, we observed that compared with control 
subjects, patients with RRMS had significantly reduced structural 
connections and increased functional connectivity of pair-wise 
subregions over long distances in the DMN. Specifically, two 
significant correlations between the FC and SC measures were found 
in the inter-subregions within the DMN. In addition, relationships 
between MS-related structural abnormalities and clinical markers 
were detected more often than were those between functional 
abnormalities and clinical markers, suggesting that the long-distance 
connection provided by the white matter is the stable foundation for 
the function of the entire network. This study, linking the increased 
inter-subregion functional coherence and decreased SC (long white 
matter damage) in pair-wise DMN subregions, provides new insights 
into the understanding of the functional reorganization in RRMS. 
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(TIF) 

References 

1. Dzicdzic T, Metz I, Dallcnga T, Konig FB, Miillcr S, ct al. (2010) Wallerian 
degeneration: a major component of early axonal pathology in multiple sclerosis. 
Brain Pathol 20: 976-985. 

2. GeurtsJJ, Barkhof F (2008) Grey matter pathology in multiple sclerosis. Lancet 
Neurol 7: 841-851. 

3. Vcrcellino M, Mascra S, Lorenzatti M, Condcllo C, Mcrola A, ct al. (2009) 
Dcmyelination, inflammation, and neurodegeneration in multiple sclerosis deep 
gray matter. J Neuropathol Exp Neurol 68: 489-502. 

4. Bruck W (2005) The pathology of multiple sclerosis is the result of focal 
inflammatory dcmyelination with axonal damage. J Neurol 252 Suppl 5: v3— 9. 

5. Bodini B, Battaglini M, De Stefano N, Khaleeli Z, Barkhof F, ct al. (201 1) T2 
lesion location really matters: a 10 year follow-up study in primary progressive 
multiple sclerosis. J Neurol Neurosurg Psychiatry 82: 72-77. 

6. Zhou F, Shiroishi M, Gong H, Zee CS (2010) Multiple sclerosis: hyperintense 
lesions in the brain on Tl-weighted MR images assessed by diffusion tensor 
imaging. J Magn Rcson Imaging 31: 789-795. 

7. Ramli N, Rahmat K, Azmi K, Chong HT (2010) The past, present and future of 
imaging in multiple sclerosis. J Clin Neurosci 17: 422^27. 

8. Rooscndaal SD, Schoonheim MM, Hulst HE, Sanz-Arigita EJ, Smith SM, et al. 
(2010) Resting slate networks change in clinically isolated syndrome. Brain 133: 
1612-1621. 



Figure S2 Example of the long white matter fibers of 
pair-wise subregions within the default-mode network 
was detected in one control subject. In this figure, the 
number indicates the subject of detected fibers by probabilistic 
tractography in the two groups. (PCC/PCUN = posterior cingu- 
late cortex/precuneus, IPL = inferior parietal lobule, mTL = me- 
dial temporal lobe, MPFC = medial prefrontal cortex, 
HC = healthy control, RRMS = relapsing remitting multiple 
sclerosis. Same abbreviated for all figure and tables). 
(TIF) 

Figure S3 An example of the relationship between 
functional connectivity and structural connectivity mea- 
sures in others threshold (0.05, 0.10, 0.15, 0.20, 0.25, 
and 0.30) of probabilistic tractography. 

(TIF) 

Table SI The union of default mode network subre- 
gions from healthy control subjects and RRMS patients 
(P<0.05, false discovery rate corrected). 

(DOC) 

Table S2 Measures of the structural connectivity (SC) 
and functional connectivity (FC) of pair-wise default 
mode subregions between two groups (mean ± standard 
deviation) . 

(DOC) 

Table S3 Comparison of the structural connectivity 
(SC) and functional connectivity (FC) of pair-wise 
default mode subregions between two groups [one- 
way repeated-measure analysis of variance, F values (P 
values)] . 
(DOC) 

Acknowledgments 

The authors thank their patients and volunteers for participating in this 
study and Dr. Wei Liao and Prof. Yufeng Zang of Hangzhou Normal 
University for their assistance. 

Author Contributions 

Conceived and designed the experiments: FZ HG. Performed the 
experiments: BW XW QC. Analyzed the data: BW FZ XW LW. 
Contributed reagents/ materials/ analysis tools: FZ BW YZ HW. Wrote the 
paper: FZ. 



9. Faivre A, Rico A, Zaaraoui W, Crcspy L, Rcuter F, ct al. (2012) Assessing brain 
connectivity at rest is clinically relevant in early multiple sclerosis. Mult Scler 18: 
1251-1258. 

10. Audoin B, Au Duong MV, Malikova I, Confort-Gouny S, Ibarrola D, ct al. 
(2006) Functional magnetic resonance imaging and cognition at the very early 
stage of MS. J Neurol Set 245: 87-91. 

11. Rocca MA, Valsasina P, Absinta M, Riccitclli G, Rodeghcr ME, ct al. (2010) 
Default-mode network dysfunction and cognitive impairment in progressive MS. 
Neurology 74: 1252-1259. 

12. Koch MA, Norris DG, Hund-Gcorgiadis M (2002) An investigation of 
functional and anatomical connectivity using magnetic resonance imaging. 
Neuroimage 16: 241-250. 

1.3. Grcicius MD, Supckar K, Mcnon V, Dougherty RF (2009) Rcsting-statc 
functional connectivity reflects structural connectivity in the default mode 
network. Ccreb Cortex 19: 72-78. 

14. van den Heuvel MP, Mandl RC, Kahn RS, Hulshofl' Pol HE (2009) 
Functionally linked resting-state networks reflect the underlying structural 
connectivity architecture of the human brain. Hum Brain Mapp 30: 3127-3141. 

15. Honey CJ, Thivierge JP, Sporns O (2010) Can structure predict function in the 
human brain? Neuroimage 52: 766-776. 

16. Schoonheim MM, Gcurts JJ, Barkhof F (2010) The limits of functional 
reorganization in multiple sclerosis. Neurology 74: 1246-1247. 



PLOS ONE | www.plosone.org 



8 



July 2014 | Volume 9 | Issue 7 | e101 198 



Altered Inter-Subregion Connectivity in the DMN 



17. Schoonheim MM, Filippi M (2012) Functional plasticity in MS Friend or foe? 
Neurology 79: 1418-1419. 

18. Bonavita S, Gallo A, Sacco R, Cortc MD, Bisccco A, ct al. (2011) Distributed 
changes in default-mode resting-state connectivity in multiple sclerosis. Mult 
Scler 17: 411-422. 

19. Hawellek DJ, Hipp JF, Lewis CM, Gorbetta M, Engcl AK (2011) Increased 
functional connectivity indicates the severity of cognitive impairment in multiple 
sclerosis. Proc Nad Acad Sci USA 108: 19066-19071. 

20. Allen EA, Erhardt EB, Damaraju E, Gruner W, Segall JM, et al. (2011) A 
baseline for the multivariate comparison of resting-state networks. Front Syst 
Neurosci 5: 2. 

2 1 . Dogonowski AM, Siebner HR, Soelberg Sorensen P, Paulson OB, Dyrby TB, et 
al. (2013) Resting-state connectivity of pre-motor cortex reflects disability in 
multiple sclerosis. Acta Neurol Scand 128: 328—335. 

22. Lowe MJ, Bcall EB, Sakaie KE, Kocnig KA, Stone L, et al. (2008) Resting state 
sensorimotor functional connectivity in multiple sclerosis inversely correlates 
with transcallosal motor pathway transverse diffusivity. Hum Brain Mapp 29: 
818-827. 

23. Loitfelder M, Filippi M, Rocca M, Valsasina P, Ropele S, et al. (2012) 
Abnormalities of resting state functional connectivity are related to sustained 
attention deficits in MS. Plos One 7: e42862. 

24. Filippi M, Agosta F, Spinelli EG, Rocca MA (2012) Imaging resting state brain 
function in multiple sclerosis. J Neurol 260: 1709-1713. 

25. Andre ws-HannaJR, ReidlerJS, Huang C, Buckncr RL (2010) Evidence for the 
default network's role in spontaneous cognition. J Neurophysiol 104: 322-335. 

26. Zhang D, Raichle ME (2010) Disease and the brain's dark energy. Nat Rev 
Neurol 6: 15-28. 

27. Cruz-Gomez AJ, Ventura-Campos N, Belenguer A, Avila C, Forn C (2013) The 
link between resting-state functional connectivity and cognition in MS patients. 
Mult Scler [Epub ahead of printj. 

28. Enzinger C, DeLuca J (2012) Large-scale neuronal network dysfunction in 
multiple sclerosis? evidence from resting-state fMRI. Neurology 79: 1416-1417. 

29. van den Heuvel MP, Kahn RS, Goni J, Sporns O (2012) High-cost, high- 
capacity backbone for global brain communication. Proceeding's of the National 
Academy of Sciences 109: 11372-11377. 

30. Taylor PA, Saad ZS (2013) FATCAT: (an efficient) functional and tractographic 
connectivity analysis toolbox. Brain Connect 3: 523—535. 

31. Khalsa S, Mayhew SD, Chechlacz M, Bagary M, Bagshaw AP (2013) The 
structural and functional connectivity of the posterior cingulate cortex: 
comparison between deterministic and probabilistic tractography for the 
investigation of structure -function relationships. Neuroimage [Epub ahead of 
printj 

32. Liao W, Zhang Z, Pan Z, Mantini D, Ding J, et al. (20 1 1) Default mode network 
abnormalities in mesial temporal lobe epilepsy: a study combining fMRI and 
DTI. Hum Brain Mapp 32: 883-895. 

33. Horn A, Ostwald D, Reisert M, Blankenburg F (2014) The structural- functional 
conneetome and the default network of the human brain. Neuroimage. [Epub 
ahead of print] 

34. Polman CH, Reingold SC, Banwcll B, Clanet M, Cohen JA, et al. (2011) 
Diagnostic criteria for multiple sclerosis: 2010 revisions to the McDonald 
criteria. Ann Neurol 69: 292-302. 

35. Hurwitz BJ (2009) The diagnosis of multiple sclerosis and the clinical subtypes. 
Ann Indian Acad Neurol 12: 226-230. 

36. Kurtzke JF (2008) Historical and clinical perspectives of the expanded disability 
status scale. Ncurocpidcmiology 31: 1-9. 



37. Shu N, Liu Y, Li K, Duan Y, Wang J, et al. (201 1) Diffusion tensor tractography 
reveals disrupted topological efficiency in white matter structural networks in 
multiple sclerosis. Cereb Cortex 21: 2565-2577. 

38. Pclletier D, Garrison K, Henry R (2004) Measurement of whole-brain atrophy 
in multiple sclerosis. J Neuroimaging 14: 11S-19S. 

39. Yan CO, Cheung B, Kelly C, Colcombe S, Craddock RC, et al. (2013) A 
comprehensive assessment of regional variation in the impact of head 
mieromovements on functional conncctomics. Neuroimage 76: 183—201. 

40. Li YO, Adali T, Calhoun VD (2007) Estimating the number of independent 
components for functional magnetic resonance imaging data. Hum Brain Mapp 
28: 1251-1266. 

41. Qi R, Xu Q, Zhang LJ, Zhong J, Zheng G, et al. (2012) Structural and 
functional abnormalities of default mode network in minimal hepatic 
encephalopathy: a study combining DTI and fMRI. Plos One 7: c41376. 

42. Jenkinson M, Bannister P, Brady M, Smith S (2002) Improved optimization for 
the robust and accurate linear registration and motion correction of brain 
images. Neuroimage 17: 825—841. 

43. Behrens TE, Berg HJ, Jbabdi S, Rushworth MF, Woolrich MW (2007) 
Probabilistic diffusion tractography with multiple fibre orientations: What can 
we gain? Neuroimage 34: 144—155. 

44. Morris DM, Embleton KV, Parker GJ (2008) Probabilistic fibre tracking: 
differentiation of connections from chance events. Neuroimage 42: 1329-1339. 

45. Zhou F, Zee CS, Gong H, Shiroishi M, LiJ (2010) Differential changes in deep 
and cortical gray matters of patients with multiple sclerosis: a quantitative 
magnetic resonance imaging study. J Comput Assist Tomogr 34: 431. 

46. Roosendaal S, Geurts J, Vrenken H, Hulst H, Cover K, et al. (2009) Regional 
DTI differences in multiple sclerosis patients. Neuroimage 44: 1397—1403. 

47. Sbardella E, Tona F, Petsas N, Pantano P (2013) DTI measurements in multiple 
sclerosis: evaluation of brain damage and clinical implications. Mult Scler Int 
2013: 671730. 

48. Ciccarelli O, Wheeler-Kingshott CA, McLean MA, Cercignani M, Wimpey K, 
et al. (2007) Spinal cord spectroscopy and diffusion-based tractography to assess 
acute disability in multiple sclerosis. Brain 130: 2220—2231. 

49. Kern KC, SarconaJ, Montag M, Giesser BS, Sicotte NL (2011) Corpus callosal 
diffusivity predicts motor impairment in relapsing-remitting multiple sclerosis: a 
TBSS and tractography study. Neuroimage 55: 1169-1 177. 

50. Gorgoraptis N, Wheeler-Kingshott CA, Jenkins TM, Altmann DR, Miller DH, 
et al. (2010) Combining tractography and cortical measures to test system- 
specific hypotheses in multiple sclerosis. Mult Scler 16: 555-565. 

51. Hu B, Ye B, Yang Y, Zhu K, Kang Z, et al. (201 1) Quantitative diffusion tensor 
deterministic and probabilistic fiber tractography in relapsing-remitting multiple 
sclerosis. EurJ Radiol 79: 101-107. 

52. Joel SE, Caffo BS, van Zijl PC, Pekar JJ (2011) On the relationship between 
seed-based and ICA-bascd measures of functional connectivity. Magn Reson 
Med 66: 644-657. 

53. Kaiser M, Hilgetag CC (2006) Nonoptimal component placement, but short 
processing paths, due to long-distance projections in neural systems. PLoS 
Comput Biol 2: c95. 

54. van Oort ES, van Gappellen van Walsum AM, Norris DG (2013) An 
investigation into the functional and structural connectivity of the Default 
Mode Network. Neuroimage 90C: 381-389. 

55. Ciccarelli O, Werring D, Wheeler-Kingshott C, Barker G, Parker G, et al. 
(2001) Investigation of MS normal- appearing brain using diffusion tensor MRI 
with clinical correlations. Neurology 56: 926-933. 

56. Filippi M, Riccitelli G, Mattioli F, Capra R, Stampatori C, et al. (2012) Multiple 
sclerosis: effects of cognitive rehabilitation on structural and functional MR 
imaging mcasurcs-an explorative study. Radiology 262: 932-940. 



PLOS ONE | www.plosone.org 



9 



July 2014 | Volume 9 | Issue 7 | e101 198 



